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Optimal engine performance using inference for non-identical finite source 
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We quantify the prior information to infer the optimal characteristics for a constrained 
thermodynamic process of maximum work extraction for a pair of non-identical finite sys¬ 
tems. The total entropy of the whole system remains conserved. The ignorance is assumed 
about the final temperature of the finite systems and then a prior distribution is assigned to 
the unknown temperatures. We derive the estimates of efficiency for this reversible model 
of heat engine with incomplete information. The estimates show good agreement with effi¬ 
ciency at optimal work for arbitrary sizes of systems, however the estimates become exact 
when one of the reservoir becomes very large in comparison to the other. 
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I. INTRODUCTION 


In this paper, we revisited the problem of maximum work extraction from the inference ap¬ 
proach [l-^. The extraction of maximum work using a pair of finite source and sink has been 
discussed earlier in literature [b-lO]. The whole set-up working as a heat engine delivers maxi¬ 
mum work for a reversible process in accordance with Maximum Work Theorem. The maximum 
achievable efficiency of such a heat engine is clearly is less than that of the Carnot efficiency 
[rjc) because of the finite size of source/sink. In recent years, the field of finite-time thermo¬ 
dynamics has become a popular area of research since it deals with the realistic constraints in 
heat engines such as finite-time of the engine operation cycle, finite reservoirs, internal friction 
etc. These constraints, in turn, lead to a less efficient heat engine than Carnot engine but is of 
practical importance. However, such engines can be optimised to deliver maximum power with 
a compromise of efficiency. 
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The inference approach is based on Bayesian approach to uncertainty 


11I414|| in which any 


uncertainty can be treated probabilistically. This probabilistic approach is associated with a 
rational de gree of belief of the observer rather than relative frequency interpretation of prob¬ 
ability |l5l . Ilo ]. The uncertain parameter is assigned with a probability distribution simply 
known as a prior to quantify the ignorance of the likely values of the uncertain parameter. 
Priors with different justification have been proposed [l2, 17-201. In our work, we propose prior 
for a constrained thermodynamic process with incomplete information of the thermodynamic 
coordinates. 

In our previous work Q] , we addressed the problem of maximum work extraction with finite 
source/sink within this inference approach. A pair of identical, similar sized hnite reservoirs 
were considered which serve as finite source and sink. The fundamental thermodynamic relation 
obeyed by the reservoirs is taken to be S' = , where k may depend on some universal con¬ 

stants and/or volume, particle number of the system. We restrict to the case 0 < wi < 1, which 
implies systems with a positive heat capacity. The optimal or maximum work extracted from 
this reservoir set-up is estimated. The efficiency at optimal work is also inferred upto second 


order as r] ^ + 'He'/‘^ 


2ll-l24i|. in near equilibrium regime {rje <C 1). A generalisation of this 


n 


^proach can be thought of by considering the non-identical systems as reservoirs [^. In paper 
[4], finite reservoirs are modelled by perfect gas systems with different constant heat capacities. 
Thus, the new information about distinct source and sink was utilized in the assignment of 
prior for the uncertain temperatures. The temperatures Ti and T 2 can be distinguished now. 
Moreover, the ranges of allowed values of Ti and T 2 are different. In this paper, we consider two 
dissimilar systems obeying a thermodynamic relation of the form S = kU‘^^ and reconsider the 
maximum work extraction process within inference approach. We derive the temperature and 
efficiency estimates for this model which show remarkable agreement with their optimal values. 

This paper is organised as follows. In section m we discuss the model for finite reservoirs. 
Section m outlines the discussion of the range for Ti and T 2 and the form of priors. It also 
comprises of the discussion of inference for special cases when one system becomes very large 
in comparison to the other. Then, we discuss the estimation procedure close to equilibrium 
analytically. In section IV A[ numerical results for arbitrary sizes of reservoirs have discussed. 
Finally in section ED we make some concluding remarks on our extended inference approach 
applied in case of non-identical systems. 
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II. MODEL 


To model the finite reservoirs, consider a pair of thermodynamic systems obeying the relation 
of the form S = where U is the internal energy of the system and cui is some known 

constant. Some well-known physical examples in this framework are the ideal Fermi gas (wi = 


1/2), the degenerate Bose gas (wi = 3/5) and the black body radiation (wi = 3/4) 


25l |. Classical 


ideal gas can also be treated as the limit, oji —>■ 0. Using the basic definition : dS/dU = 1/T, 
we get: U = (wikT)^^‘^. Alternately, we can write : S = {coiT)‘^, where oj = wi/(l — wi). 

Since the thermodynamic relation obeyed by the two systems remains the same, the two 
may be non-identical only if they differ in their volumes, number/nature of particles etc. Thus, 
it is the constant of proportionality, k, which is different for the two systems. Let T+ and 
T_ (< T+) be the initial temperatures of two systems with ni and K 2 as the proportionality 
constants respectively. 

To perform inference, examine an arbitrary intermediate stage of the process when the 
temperatures of the two systems are T[ and T 2 respectively. The work extracted from the 
engine is W = — AU which can be written as: 


W = (k2CJiy^^ 




( 1 ) 


where 7 = K 2 /^i- For convenience, we define 7 ^+“^ = cr, 6 = T_/T+, Ti = r)/T+ and T 2 = 
T 2 /T-. Thus, work can be rewritten as : 


W = (k2UJiT+)^+‘^ - Ti^+‘^) + - T 2 ^+‘^)] . 


(2) 


The constraint of entropy conservation AS = 0 gives Si + S 2 = S-f. -h S-, where 5+ (S-) and 
(S 2 ) are the values of entropy in initial and final state of system 1 ( 2 ) respectively. Then, 
we can write: 


Ti = ll + a(0‘^-T2‘^)]^ , 


( 3 ) 


or equivalently. 


T2 = + 


( 4 ) 


By making use of above equations, work can be written as a function of one variable (say T 2 ) 
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only: 


W{T2) 


(fe2WlT+)l+- 



(1 + ^( 0 - 


72 '^))^) + [e^+^ -T2^+‘^) 


(5) 


Similarly, work can be written as a fnnction of Ti also. 

The optimal work can be extracted from the engine when the two systems reach a common 
temperature. The common temperature, Tc , is given from Eq. ([3]) with Ti = T 2 = as: 

The efficiency of the engine is given as r/ = 1—Qout/Qin, where Qin = 
is the heat absorbed from the source and Qout = — 0 ^+'^) is the waste heat 

rejected to the sink. Thus, for any arbitrary value of 7 , efficiency at any intermediate stage of 
the process can be given as: 


rj^ = 1 — a 


(7^21+'^ - 
(1 - 


(7) 


For efficiency at optimal work {rj*), we substitute Ti = T 2 = Tc in above equation to obtain: 


r]* = l-a 


(i-rci+‘") 


( 8 ) 


Let us discuss the efficiency in the limiting cases: 

(a) In the limit 7 —?> 0 i.e. when heat source is very large as compared to heat sink, temperature 
of source remains constant at = T+ or Ti = 1 while temperature of sink approaches this 
value for optimal work extraction. We write efficiency as a function of T 2 as 


uj \ fT2^+^ 

l + u) V T2^ - 0^ J ■ 


(9) 


Efficiency at optimal work in this limit is given by substituting r 2 = 1 in above equation as: 


* 1 
% = 1 - 


UJ 


1 + w 


1 - 
1-9^^ 


( 10 ) 


(b) In the limit 7 —00 i.e. when heat sink is very large in comparison to heat source, the sink 
stays at temperature (T_) and source approaches this value for optimal work extraction. The 


4 











efficiency can be expressed in terms of Ti as: 


r/oo = 1 - 


l + wN { 1-Ti 


UJ 


1-Ti 


l+o; 


For the optimal process, substitute Ti = 0 in above expression to obtain: 


= 1 - ^ 


l + a;\ / 1- 


UJ 


1 - / ■ 


III. ASSIGNMENT OF PRIOR 


( 11 ) 


( 12 ) 


The inference is performed by assigning the prior probability distributions for the uncertain 
parameters Ti and T2, since we assume ignorance of the actual values of Ti and T2 or the 
extent to which the process has proceeded. Thus, there are two observers, 1 and 2 , for Ti 
and T2 respectively. Let us hrst summarise the prior information we possess before making the 
inference: 

(i) There exists a one-to-one relation between Ti and T2 given by Eq. ([ 3 ]) or Eq. (| 3 |) which 
suggests that probability of Ti to lie in small range \Ti,Ti + dTi] is same as the probability of 
T2 to lie in [r2, r2 -|- dT2\, so we can write : 


Pi{Ti)dTi = P2{T2)dT2, 


( 13 ) 


where Pi and P2 are the normalised prior distribution functions for Ti and T2. 

(ii) The set-np works like a heat engine and thus W > 0 . 

With identical systems, we had an additional assumption of the same form of normalised 
prior distribution, P, for Ti and T2 and thus prior distribution is invariant under the change of 
parameter. Here, with non-identical systems, still we assume the same functional form / of the 
prior distribntions for Ti and T2. However, since now the two systems are not identical, this 
information has to be incorporated while assigning the prior. So, this information is incorporated 
in the allowed values of the range of Ti and T2. Now, the allowed range for Ti and T2 is not 


[ 6 , 1 ]. It will be different for both the parameters for different values of 7 (7^ 1 ). Thus 


say, Ti 


ranges in [Tm, 1 ] and T 2 ranges in [6,Tm\ respectively satisfying the constraint VF > 0 [^. This 
can be shown graphically in Figured] for the two cases with 7 < 1 and 7 > 1 . For 7 < 1 (larger 
source in comparison to sink), the range of allowed values of Ti is narrower than the range for 
T2. In the limiting case of 7 —>■ 0 (infinite source and finite sink), [Tm,!] shrinks to a point 
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T2{T,) 71 ( 72 ) 


FIG. 1: Work as a function of Ti (Dashed Curve) and T 2 (Solid Curve) for ui = 0.1, 6 = 0.4 


Ti = 1 which is expected for an infinite source as now the temperature of the source stays at 
r+ = i. Similarly, for 7 > 1 (larger sink in comparison to source), the range for T 2 shrinks in 
comparison to the range for Ti, and [0,Tm\ shrinks to a point T 2 = 9 for 7 —>■ 00 (infinite sink 
and finite source). This information on the range of uncertain parameters will be incorporated 
to determine the normalisation constants for prior distributions and thus we can write Pi and 
P 2 as ^\: 

PiiTi) 

P2{T2) 

where the form of / is determined from the constraint condition dS = 0, which can be written 
as: 


^ /(Ti) 

/^^/(ri)dri’ 

^ fm 

f{T2)dT2 


dSi + dS2 = 0 , 


which can be further written, using additivity of energy, as: 


dSi 


V 


m] 

dPj 


dTi + 
V 




dT2 = 0 . 

V 


Using {^)v = ^ and (fy)y = C{T) in above equation, we get: 


m ^ C2IT2 
dT2 Ci/Ti^ 


(16) 


(17) 


(18) 
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Combining Ea. ffTSll and Ea. (fT3]l . suggests the following form of prior 


fl: 


P{T^) = 


Ci{Ti)lTi 


N 


(19) 


where i = 1, 2 and N = fg Ci{Ti)/Ti dTi is the normalisation constant. With our model of the 
reservoirs {Si oc the functional form of the prior distribution can be written as: 


PiT,) oc Ti 


(jJ—1 


( 20 ) 


IV. ESTIMATION OF TEMPERATURE 


The expected value of a temperature is: 


r-T, 


Ti = 


iTi, 


TiP{Ti)dTi, 


( 21 ) 


where i = 1,2. Taking into account the respective ranges of allowed values of Ti and T 2 , 
identified above, we obtain: 


Ti = 


UJ 


1 -\- UJ 


l-Tr. 


l+o; 


1-T„ 


( 22 ) 


and 


T2 = 


OJ 


^ _/3a; 


l+uj \ Tm^ -9' 


(23) 


To determine Tm or Tm-, we solve Eq. ([2|) by setting W{Ti) = 0 or W{T 2 ) = 0 respectively. In 
general, Eq. (ED has to be solved numerically for arbitrary values of coi. For ideal Fermi gas 
(wi = 1/2), it can be solved analytically [2^ and thus we get: 


Ti € 
T2 G 


1-72 + 2720 ■ 

,2 ’ ^ 


9, 


1 + 7' 
2-0 + 720 
1 + 72 


(24) 

(25) 


Due to Eq. we can write one-to-one relation between Tm and Tm as: 


1 - = a(rM" - 0‘"). 


(26) 
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Using above equation in W{Tm) = 0 , we obtain: 


1 - (27) 

From Eqs. ([22]), (f23]l . (1261) and (l27|l . we can write: 

Ti = Ta. (28) 

However, firstly we will solve the Eq. m for the limiting cases when one of the systems become 
very large in comparison to the other system. 


A. Infinite source and finite sink 


This case corresponds to the limit 7 —>• 0. Here the only uncertain parameter is T 2 as 
temperature of source stays at r+ = 1 while the temperature of sink approaches r+ at optimal 
work extraction. To discuss this limit, we set Eq. ([5|) as W{T 2 ) = 0 and obtain: 


rp^-\-U} _ _ — 

(J 


1 - (1 + a{9^ - T2^)) 


l + CJ 


(29) 


Taking the limit 7 ^ 0 , the above equation gets simplified to : 


a; {T 2 ^+^ - 9^+^) = (1 + w) (Ts"^ - 9 ^), 


(30) 


whose trivial solution is T 2 = 0. The other solution is Tm so we write: 

w - 9^+^) = (1 + w) {Tm^ - r). (31) 

Consistency between Eqs. (f23|) and (f3T]l demands that we must have: 


T2 = 1. 


(32) 


Thus expected sink temperature exactly matches with temperature of heat source for optimal 
process. The efficiency is estimated by replacing T 2 in Eq. (I9|) by Eq.([32|) and estimate for 
efficiency is same as Eq. m- Hence, inference approach reproduces the optimal behaviour 
exactly in the limit 7 —>■ 0 . 




B. Finite source and infinite sink 


Consider the case of infinite sink in comparison to source (7 —>■ 00 ). Here, the sink stays at 
temperature T-{= 9) and the temperature of source approaches T_ for optimal work extraction. 
Hence Ti is the only uncertain parameter for this limiting case. The range for Ti is determined 
by using ([1|) in ([2]) and then setting W{Ti) = 0, we get: 


= a 


(O' 




1 + CJ 


-9 


l+CJ 


(33) 


In the limit j —>■ 00 , the above equation gets simplifies to: 


w(l-ri^+‘^) =0(l + w)(l-Ti‘^). (34) 

The trivial root of above equation is Ti = 1 and other root (T^) satisfies: 

cu (1 - = 9(1 + w) (1 - T^-). (35) 

From Eqs. (f2^ and (f35]l . we obtain: 


Ti = 9. (36) 

It is clear from the Eq. (|36|) that the average temperature of the source exactly matches with the 
temperature of the infinite sink which happens in case of maximum work extraction. Eurther, 
efficiency at optimal work ( 7 ^) is also inferred exactly due to Eqs. (fTT]l and (f36]l . Thus, we 
are able to infer exactly the optimal behaviour of the system with infinite sink and finite source 
also. 


C. Near-equilibrium estimation 

In this section, we approximate the values of and Tm when 9 is close to unity. Eor this, 
consider the case 0 < 7 < 1. Let us examine the case of observer 2 . Since close to equilibrium, 
Tm is also close to unity so we can introduce a small parameter e > 0 such that Tm = 0 (1 + e). 
Rewriting the Eq. dS]) as W(Tm) = 0 : 

1 + cj(9^+‘^[l - (1 + 6)^+^^] = (1 + a9^[l - (1 + e)^])^ . (37) 
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Making series expansion in e and keeping terms only upto second order, we get a quadratic 
equation in e as: 


(1 - cu ) - w(l - 9 )+ 2 ] 

- 3 [ o -- a ;( l - 6 l ) + l]e + 6 ( l - 6 l ) = 0 . ( 38 ) 


For instance, if we take limit a; —>■ 0 in above equation, we reproduce the case for perfect gas as: 


(7 + 1)(7 + 2 ) e ^ - 3(7 + l)e + 6(1 - 9 ) = 0 , 


whose acceptable solution 


Q 


is approximated upto second order in rjc as: 


e 


2 

1 + 7 


7c + 


i ( 2 ± 7 _ , 
3 ( 1 + 7 )= 


Similarly, Eq. (l38)) can be solved for e and the solution can be approximated as: 


(39) 


(40) 


e 


^ 2 [ A-uj + a{uj + 2 )] 2 

1 + u ^" 3 (l + a )2 


(41) 


Then Tm is determined, which in turn determines due to Eq. ([5|). 

Suppose Ti (T 2 ) are the estimates for Ti (r 2 ) by the observer 2 (1) by making use of (l3|) and 
©• Let us examine the near-equilibrium expansion of the estimates of temperature as well as 
Tc as: 


„ a cj(l — cj) 9 q, 

^2 = 1-r—^7 c- 97-i , _79 7e+0[7c]l 


l + u '^- 3 ( 1 + a )^ 


T 2 = 1- 


£7 (1 + 3 cr)(l - oj ) 2 , 31 


1 + cr 


7 c 


6(1 + it)2 


- 7 c + 0 [ 7 j . 


Ti = 1- 


a a{3 + a){l - w) 2 , 3i 


1 + cr 


7 c 


6(1 + ( t )2 


- 7 c + 0[Vc 


Tr = I- 


(7 


1 + 2(1 + £ 7)2 


£ 7 ( 1 - 0 ;) 2 

'/c 


£ 7(2 + £7 + o ;(£ 7 - 1 ))( 1 - o ;) 3 4 

6(1 + £ 7)3 ^ 


(42) 


(43) 


(44) 


(45) 


We have skipped the lengthy expressions of third order terms in series expansion of temperature 
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estimates. Let us define the estimated value of the temperature of one reservoir as the weighted 
mean of the estimates by two observers. It will be given as 


Tim — ^Ti + BTi, 


(46) 


where i = 1,2 and A, B are the weights satisfying the condition A + B = 1. In the above case, 
we have seen that both estimates, T 2 (by observer 2 ) and T 2 (by observer 1 ), match with 
only upto first order. If we choose weights A and B so as to obtain matching beyond first order, 
then the weights are calculated as: 


A 

B 


1 

1 + fj’ 

a 

1 + (T 


(47) 

(48) 


This weighted mean estimated temperature {T 2 ^m) of the sink shows remarkable agreement with 
Tc upto third order in rjc close to equilibrium. Further, for 7 —)■ 0, T 2 ^m becomes exactly equal 
to r 2 showing that estimation is done only by observer 2 as the other reservoir corresponding 
to observer 1 (source) becomes infinite in comparison and hence, its temperature stays constant 
at T+ and no uncertainty exists in its value as the process proceeds. In other limit of 7 —>■ 00 , 
Ti becomes exactly equal to Tc- 


V. ESTIMATION OF EFFICIENCY 

Efficiency is estimated by replacing r 2 in Eqs. m and dSI) to obtain efficiency estimate by 
observer 2 {if 2 ). Similarly, denote the efficiency estimate by observer 1 as ( 71 ). Expanding the 
estimates of efficiency close to equilibrium and make a comparison with the efficiency at optimal 
work as follows: 


Vi 

m 

* 

^7 


,,, 2(1 + a) + u(a- 2) , , 

2 ^ 12(1 +a) 

(49) 

nc (i + +)++(2<7-i), 3 

2 ^ 12(1 +a) 

(50) 

,k (1 + 2a) + +(a-l), 3 

2 12(1 + a) 

(51) 


It is clear from the above expressions that estimates of efficiency either by observer 1 (fji) or 
by observer 2 ( 72 ) matches with efficiency at optimal work (r/*) only upto first order in ijc- 


11 







However, we define mean efficiency {rim) and compare it with rj* as: 


fim = Afji + H? 72 , (52) 

where A and B are the weights assumed for the estimation done by observer 1 and observer 
2 respectively satisfying A + B = 1, similar to the case of temperature estimation. Estimated 
mean efficiency (? 7 m) matches with efficiency at optimal work upto second order. For the limiting 
cases of 7 —)■ 0 and 7 —>■ oo, 72 and fji give the exact estimates for efficiency at optimal work 
respectively. 


A. Numerical results for arbitrary 7 


Eq. ([5]) has to be solved to determine the roots, one trivial root is 0 while the other root, 
Tm, can be determined numerically for given values of wi, 7 , and 9. Then, we obtain numerical 
estimates of efficiency by observer 1 ( 71 ) and observer 2 ( 72 ) for arbitrary values of 7 . Figure 
[ 2 ] shows the comparative plots of efficiency for different 7 ’s. 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

G 6 


FIG. 2: Efficiency as a function of 6 for cji = 0.1. The dashed curve is for efficiency at 
optimal work, uppermost curve (blue online) isforpi, lowermost curve (purple online) is for 
72 , and middle curve (red online) is for mean efficiency (fjm)- 


From the numerical plots, it has become clear that in the 0 < 7 < 1 regime, estimates made 
by observer 2 give better results as compared to the observer 1 and vice-versa in the case of 

7 > 1 . 
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VI. CONCLUSION 


Thus, we have extended our previous approach of inference studied in where two identical 
finite systems were taken as source and sink. Taking non-identical systems, we can distinguish or 
label the two systems acting as source and sink. The range for Ti and T 2 is also different because 
of dissimilar systems. While generalising this approach, we have observed that estimates match 
exactly with their optimal values when one of the reservoir becomes very large as compared to 
the other. For arbitrary values of 7 also, numerical calculations have been performed. These 
calculations show that information incorporated in the prior distribution reproduce the optimal 
behaviour of the system, however, now the efficiency estimates made by two observers are not 
symmetrically distributed about the efficiency at optimal work unlike in the case with similar 
reservoirs (7 = 1). Instead, estimates made by one observer he closer to the optimal value as 
compared to the other depending upon the value of 7. Near-equilibrium, it has been observed 
that universality, r/^/8, in efficiency does not hold and becomes system dependent. However, 
in this case, efficiency at optimal work can be reproduced upto second order by defining mean 
efficiency with non-identical weights for the efficiency estimates by the two observers. Thus, 
with non-identical systems also, we quantify the prior information and use it to estimate the 
optimal performance in constrained thermodynamic process 
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